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Abstract 

Using gauge theory /gravity duality we study sound wave propagation in strongly cou- 
pled non-conformal gauge theory plasma. We compute the speed of sound and the 
bulk viscosity of M = 2* supersymmetric SU (N c ) Yang-Mills plasma at a temperature 
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1 Introduction and summary 



The conjectured duality between gauge theory and string theory [1, 2, 3, 4] is useful for 
insights into the dynamics of gauge theories. For strongly coupled finite-temperature 
gauge theories in particular, the duality provides an effective description in terms 
of a supergravity background involving black holes or black branes. A prescription 
for computing the Minkowski space correlation functions in the gauge/gravity corre- 
spondence [5, 6] makes it possible to study the real-time near-equilibrium processes 
(e.g. diffusion and sound propagation) in thermal gauge theories at strong coupling 
[7, 8, 9, 10, 11, 12, 13, 14, 15]. Computation of transport coefficients in gauge theo- 
ries whose gravity or string theory duals are currently known may shed light on the 
true values of those coefficients in QCD, and be of interest for hydrodynamic mod- 
els used in theoretical interpretation of elliptic flows measured in heavy ion collision 
experiments at RHIC [16, 17, 18]. This expectation is supported by an intriguing ob- 
servation [13, 15, 19, 20] that the ratio of the shear viscosity to the entropy density in 
a strongly coupled gauge theory plasma in the regime described by a gravity dual is 
universal for all such theories and equal to 1/Att. (Finite 't Hooft coupling corrections 
to the shear viscosity of M = 4 supersymmetric SU (N c ) gauge theory plasma in the 
limit of infinite N c were computed in [14].) 

Previously, it was shown [10] that the dual supergravity computations reproduce 
the expected dispersion relation for sound waves in the strongly coupled M = 4 super- 
symmetric Yang-Mills (SYM) theory plasma 

uj(q)=v s q-t^q 2 + 0(q 3 ), (1.1) 



(1.2) 



where 

WJ 

is the speed of sound, P is pressure and £ is the volume energy density. The attenuation 
constant T depends on shear and bulk viscosities rj and (, 

1 r C+U)- (1-3) 



£ + P V 3 

For homogeneous systems with zero chemical potential £ + P = s T, where s is the 
volume entropy density. Conformal symmetry of the M = 4 gauge theory ensures that 

v. = ^=, C = 0. (1.4) 
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Indeed, in conformal theories the trace of the stress-energy tensor vanishes, implying 
the relation between S and P of the form S = 3P, from which the speed of sound 
follows. 

Non-conformal gauge theories 1 , and QCD in particular, are expected to have non- 
vanishing bulk viscosity and the speed of sound different from the one given in Eq. (1.4). 
Lattice QCD results for the equation of state 8 = P(S) suggest that v s ~ l/VS ~ 
0.577 for T ~ 2T C , where T c m 173 MeV is the temperature of the deconfining phase 
transition in QCD. When T — > T c from above, the speed of sound decreases rather 
sharply (see Fig. 11 in [21], [22], [23], and references therein). To the best of our 
knowledge, no results were previously available for the bulk viscosity of non-conformal 
four-dimensional gauge theories 2 . 

In this paper we take a next step toward understanding transport phenomena in 
four-dimensional gauge theories at strong coupling. Specifically, we study sound wave 
propagation in mass-deformed Af = 4 SU (N c ) gauge theory plasma. In the language 
of four- dimensional Af = 1 supersymmetry, Af = 4 gauge theory contains a vector 
multiplet V and three chiral multiplets $, Q, Q, all in the adjoint representation of the 
gauge group. Consider the mass deformation of the Af = 4 theory, where all bosonic 
components of the chiral multiples Q, Q receive the same mass rrib, and all their fermion 
components receive the same mass rrif. Generically, rrif ^ rrib, and the supersymmetry 
is completely broken. When rrib = rrif = m (and at zero temperature), the mass- 
deformed theory enjoys the enhanced Af = 2 supersymmetry with V, $ forming an 
Af = 2 vector multiplet and Q, Q forming a hypermultiplet. In this case, besides 
the usual gauge- invariant kinetic terms for these fields 3 , the theory has additional 
interaction and hypermultiplet mass terms given by the superpotential 

W = ^ Tr (lQ,Q]$) + -£L (ttQ 2 + TrQ 2 ) . (1.5) 
9ym v ' 9ym v ' 

1 One should add at once that throughout the paper terms "conformal" and "non-conformal" refer 
to the corresponding property of a theory at zero temperature. Conformal invariance of Af = 4 SYM 
is obviously broken at finite temperature, but we still refer to it as "conformal theory" meaning that 
the only scale in the theory is the temperature itself. 

2 Perturbative results for the bulk viscosity and other kinetic coefficients in thermal quantum field 
theory of a scalar field were reported in [24]. Bulk viscosity of Little String Theory at Hagedorn 
temperature was recently computed in [25] using string (gravity) dual description. The speed of 
sound of a four-dimensional (non-conformal) cascading gauge theory was computed in [26] , also from 
the dual gravitational description. 

3 The classical Kahler potential is normalized according to 2/ ' g\ M Tr[$$ + QQ + QQ]. 
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The mass deformation of M = 4 Yang-Mills theory described above is known as the 
M = 2* gauge theory. This theory has an exact nonperturbative solution found in 
[27]. Moreover, in the regime of a large 't Hooft coupling, g\ M N c ^> 1, the theory has 
an explicit supergravity dual description known as the Pilch- Warner (PW) flow [28]. 
Thus this model provides an explicit example of the gauge theory/gravity duality where 
some aspects of the non-conformal dynamics at strong coupling can be quantitatively 
understood on both gauge theory and gravity sides, and compared. (The results were 
found to agree [29, 30].) Given that the M = 2* gauge theory is non-conformal, and 
at strong coupling has a well-understood dual supergravity representation, it appears 
that its finite-temperature version should be a good laboratory for studying transport 
coefficients in non-conformal gauge theories. 

The supergravity dual to finite-temperature strongly coupled M = 2* gauge theory 
was considered in [31]. It was shown that the singularity-free nonextremal deformation 
of the full ten-dimensional supergravity background of Pilch and Warner [28] allows 
for a consistent Kaluza-Klein reduction to five dimensions. These nonextremal geome- 
tries are characterized by three independent parameters, i.e., the temperature and the 
coefficients of the non-normalizable modes of two scalar fields a and x ( see Section 2 
for more details). Using the standard gauge/gravity correspondence, from the asymp- 
totic behavior of the scalars the coefficient of the non-normalizable mode of a-scalar 
was identified with the coefficient of the mass dimension-two operator in the M = 4 
supersymmetric Yang-Mills theory, i.e., the bosonic mass term m^; the corresponding 
coefficient of the x-scalar is identified with the fermionic mass term rrif. The supersym- 
metric Pilch- Warner flow [28] constraints those coefficients so that mj = rrif. However, 
as emphasized in [31], on the supergravity side mj and rrif are independent, and for 
mj ^ rrif correspond to the mass deformation of M = 4 which completely breaks super- 
symmetry. In the high temperature limits m^/T <C 1 (with rrif rrib) or mj/T < 1 
(with rrii, <C iTif), the background black brane geometry was constructed analytically. 
(For generic values of rrib, rrif, the background supergravity fields satisfy a certain sys- 
tem of coupled nonlinear ODEs [31], which in principle can be solved numerically.) 
The holographic renormalization of this theory was discussed in [12], where it was also 
explicitly checked that for arbitrary values of rrib, rrif in the limit of large N c and large 
't Hooft coupling the ratio of the shear viscosity to the entropy density in the theory 
equals 1/4-7T, in agreement with universality [13, 15, 19, 20]. 

Our goal in this paper is to compute the parameters of the dispersion relation (1.1) 
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for the M = 2* plasma. Since the ratio rj/s is known, this will allow us to determine 
the speed of sound and the ratio of bulk to shear viscosity in the M = 2* theory. 
We will be able to do it only in the high-temperature regime where the metric of the 
gravity dual is known analytically, leaving the investigation of the full parameter space 
to future work. The dispersion relation (1.1) appears as a pole of the thermal two- 
point function of certain components of the stress-energy tensor in the hydrodynamic 
approximation, i.e. in the regime where energy and momentum are small in comparison 
with the inverse thermal wavelength (cu/T <C 1, q/T <^ 1). Equivalently, Eq. (1.1) can 
be computed as the lowest (hydrodynamic) quasinormal frequency of a certain class of 
gravitational perturbations of the dual supergravity background [32]. 

Though the general framework for studying sound wave propagation in strongly 
coupled gauge theory plasma from the supergravity perspective is known [10, 32], 
the application of this procedure to a non-conformal gauge theory is technically quite 
challenging. The main difficulty stems from the fact that unlike the (dual) shear mode 
graviton fluctuations, the (dual) sound wave graviton fluctuations do not decouple from 
supergravity matter fluctuations. The interactions between various fluctuations and 
their background coupling appear to be gauge-theory specific. As a result, we do not 
expect speed of sound or bulk viscosity to exhibit any universality property similar to 
the one displayed by the shear viscosity in the supergravity approximation. In fact, we 
find that the speed of sound and its attenuation do depend on the mass parameters of 
the Af = 2* gauge theory. 

Let us summarize our results. For the speed of sound and the ratio of the shear to 
bulk viscosity we find, respectively, 



where 0j 0.9672, 01 8.001, and the ellipses denote higher order terms in rrif/T and 
m b /T. From the dependence (1.6), (1.7) it follows that at least in the high temperature 
regime the ratio of bulk viscosity to shear viscosity is proportional to the deviation of 
the speed of sound squared from its value in conformal theory, 




(1.6) 




(1.7) 




(1.8) 
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where 4 re = 3tt/3J/ 2 « 4.558 for m b = 0, and re = tt 2 /3[/16 « 4.935 for m f = 0. (Note 
that the result (1.8) appears to disagree with the estimates ( ~ rj {v 2 s — 1/3) 2 [33, 34], 
later criticized in [24].) 

The paper is organized as follows. In Section 2 we review the non-extremal (finite- 
temperature) generalization of the Pilch- Warner flow. We discuss holographic renor- 
malization and thermodynamics of M = 2* gauge theory, and determine its equation 
of state to leading order in rrib/T C 1, m//T< 1 in Section 3. The equation of state 
determines the speed of sound (1.6). Computation of the sound attenuation constant 
requires evaluation of the two-point correlation function of the stress-energy tensor, 
which is a subject of Section 4. The computation follows the general strategy of [32], 
albeit with some technical novelties. We keep the parameters rrif/T and m b /T arbi- 
trary until the very end, when we substitute into equations the analytic solution for the 
metric valid in the high-temperature regime only. Alas, even in that regime we have 
to resort to numerical methods when computing the bulk viscosity. Having found the 
sound wave pole of the stress-energy tensor correlator, we confirm the result (1.6) for 
the speed of sound, and obtain the ratio of shear to bulk viscosity (1.7) for the strongly 
coupled M = 2* gauge theory in the high temperature regime. Our conclusions are 
presented in Section 6. Details of the holographic renormalization are discussed in 
Appendix A. Some technical details appear in Appendixes B-F. 

2 Non-extremal J\f = 2* geometry 

The supergravity background dual to a finite temperature M = 2* gauge theory [31] is 

a deformation of the original AdS$ x S 5 geometry induced by a pair of scalars a and \ 

of the five-dimensional gauge supergravity. (At zero temperature, such a deformation 

was constructed by Pilch and Warner [28].) According to the general scenario of a 

holographic RG flow [35], [36], the asymptotic boundary behavior of the supergravity 

scalars is related to the bosonic and fermionic mass parameters of the relevant operators 

inducing the RG flow in the boundary gauge theory. The action of the five-dimensional 

4 In general k — k(X) is a function of the ratio A = mj/m/ which we were able to compute in two 
limits A — > and A — > oo. Assuming k(A) to be a smooth monotonic function, we find that it varies 
by ~ 8% over the whole range A <G [0, +oo). Additionally, for both finite (and small) rrib/T, mf/T wc 
verified scaling (1.8) by explicit fits. 
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gauged supergravity is 



= / 

J A 



Ms 
1 



4ttG 



where the potential 5 



5 Jm 



V= — 
16 



deV=9 [\R - 3(9a) 2 - {d X f - V] 



(2.1) 



1 t dW 
3 \ da 



+ 



dW 
dx 



is a function of a, x determined by the superpotential 



W = 



-2a 



1 



e 4a cosh(2x) • 



The five-dimensional Newton's constant is 



47T 



2 5 voiss N 2 
The action (2.1) yields the Einstein equations 



= Ud^ad^a + AdpxdvX + ^g^V , 



as well as the equations for the scalars 

- 

6 da ' 



IdV 

2 dx 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



To construct a finite-temperature version of the Pilch- Warner flow, one chooses an 
ansatz for the metric respecting rotational but not the Lorentzian invariance 6 



dsl = -c\(r) dt 2 + c\(r) (dx\ + dx\ + dx\) + dr 2 
The equations of motion for the background become 



(2.7) 



a 



X " + X' (lnc lC 3)'-i^ = 0, 
c'; + c[ (\ncl)' + ^c 1 V = 0, 
4 + c' 2 (ln Cl c2)' + ^c 2 P = 0, 



(2.f 



5 We set the five-dimensional gauged supergravity coupling to one. This corresponds to setting the 
radius L of the five-dimensional sphere in the undeformcd metric to 2. 
6 The full ten-dimensional metric is given by Eq. (4.12) in [31]. 



7 



where the prime denotes the derivative with respect to the radial coordinate r. In 
addition, there is a first-order constraint 



(«0 2 + \ (x'f ~\V- J(mc 2 )'(m Cl c 2 )' = . 



(2.9) 



x e [o, l] . 



(2.10) 



3 v/w 3 2 
A convenient choice of the radial coordinate is 

x[r) = — 
c 2 

With the new coordinate, the black brane's horizon is at x = while the boundary of 
the asymptotically AdS$ space-time is at x — 1. The background equations of motion 
(2.8) become 



c' 2 ' + 4c 2 (a') 2 --4--(c 2 ) 2 + |c 2 ( X 'f 
x c 2 3 



0, 



a" H — a' — 

x 



12 Vclx 



6(a') 2 c 2 2 x + 2(xYc 2 2 x - 3c' 2 c 2 - Q{c' 2 ) 2 x 



x" + L x' 

x 



AVclx 



6(a') 2 c 2 2 x + 2(x') 2 c 2 x - 3c' 2 c 2 - 6(c' 2 ) 2 x 



^ = 

da 

ox 



(2.11) 



where the prime now denotes the derivative with respect to x. A physical RG flow 
should correspond to the background geometry with a regular horizon. To ensure 
regularity, it is necessary to impose the following boundary conditions at the horizon, 



x 







+ 



<x(x), x(x), c 2 (x) 



Si, s 2 , s 3 



(2.12) 



where Si are constants. In addition, the condition 63 > guarantees the absence of a 
naked singularity in the bulk. 

The boundary conditions at x — 1 are determined from the requirement that the 
solution should approach the AdS^ geometry as x — > 1_: 

x^l_: |a(x), x (x), c 2 (x)|-^|o,0,oc(l- a ; 2 )- 1 / 4 |. (2.13) 

The three supergravity parameters Si uniquely determine a non-singular RG flow in 
the dual gauge theory. As we review shortly, they are unambiguously related to the 
three physical parameters in the gauge theory: the temperature T, and the bosonic 
and fermionic masses mj, mj of the M = 2* hypermultiplet components. 

General analytical solution of the system (2.11) with the boundary conditions 
(2.12), (2.13) is unknown 7 . However, it is possible to find an analytical solution in 
the regime of high temperatures. 

7 One can study the system numerically, see Fig. 2 of [31]. 
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2.1 The high temperature Pilch- Warner flow 

Differential equations (2.11) describing finite temperature PW renormalization group 
flow admit a perturbative analytical solution at high temperature [31]. The appropriate 
expansion parameters are 

5 l0 c(^) 2 «l, 5 2 oc^«l. (2.14) 

Introducing a function A(x) by 

c 2 = e A , (2.15) 
to leading nontrivial order in Si, 82 we have [31] 



(2.16) 



where 



A(x) 


= 111*3-4 


a(x) 


=61 ai(x) , 


X (x) 


=$2 X2(x) , 




Oil = (1 




X2 = (l 


= -A 


f x zdz 



71 + / dy 



daA 2 (1 - y 2 



Jo (! - z 2 ) 2 \ Jo \9y J y 



(2.17) 
(2.18) 

(2.19) 



3 Jo (1 - z 2 ) 2 y J \dy J y 
The constants 7$ were fine-tuned to satisfy the boundary conditions [31]: 

8-7T 2 8 - 37T . 

* = = (2 ' 20) 

The parameters Si are related to the parameters m^, nif, T of the dual gauge theory 



via 



1 /m 6 \ 2 

: _ [ r (l)] 2 m f (2 21) 

2nT =sJl+ 1 4s 2 1 + ±S 2 2 
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Given the solution (2.16), C2 is found from Eq. (2.15), and C\ = XC2- The transition to 
the original radial variable r can be made by using the constraint equation (2.9). At 
ultra high temperatures Si — > 0, S2 — > 0, the conformal symmetry in the gauge theory 
is restored, and one recovers the usual near-extremal black three-brane metric 

dsl = (2vrT) 2 (l - x 2 )- 1 / 2 (-x 2 dt 2 + dxj + dx\ + dxf) + d f , (2.22) 

(1 x ) 

describing a gravity dual to a finite-temperature Af — 4 SYM in flat space. 



3 J\f = 2* SYM equation of state and the speed of sound 

To determine the equation of state of Af = 2* gauge theory, one needs to compute 
energy and pressure, given by the corresponding one-point functions of the stress- 
energy tensor. In computing the one-point functions, one has to deal with divergences 
at the boundary of the asymptotically AdS space which are related to UV divergences 
in the gauge theory. The method addressing those issues is known as the holographic 
renormalization [37, 38, 39, 40]. Some details of the holographic renormalization for the 
Af = 2* gauge theory are given in Appendix A. The method works for arbitrary values 
of m,b/T, m,f/T, once the solution to Eqs. (2.8) is known. In the high-temperature 
limit, energy density and pressure can be computed explicitly 



£ =-it 2 N 2 T 4 
8 

P =-n 2 N 2 T 4 
8 



64 



8 



1 + - (ln(TrT) - 1) 5f - — 6. 



71 

192 



3tt 



71 z 



ln(7rT) 5j --S. 



(3.1) 



where Si, S2 are given by Eqs. (2.21). One can independently compute the entropy 
density of the non-extremal Pilch- Warner geometry [31] , 



-w 2 N 2 T 3 



o -1 -2 
7T 7T 



(3.2) 



and verify that the thermodynamic relation, 



£-Ts = -P, 



(3.3) 



is satisfied. Alternatively, the free energy can be computed as a renormalized Euclidean 
action. Then it can be shown [12] that the free energy density T = —P obeys T = £ — 
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Ts for arbitrary mass deformation parameters m&/T, m,f/T. Finally, using Eq. (2.21) 
it can be verified that 

dS = Tds . (3.4) 

These checks demonstrate that the Af = 2* thermodynamics is unambiguously and 
correctly determined from gravity. 

We can now evaluate leading correction to the speed of sound in Af = 2* gauge 
theory plasma at temperatures much larger than the conformal symmetry breaking 
scales rrif, m^. Using Eqs. (3.1), (2.21) we find 

2 OP dP fdSY 1 1/ 64 r2 8 r2 \ 

where ellipses denote higher order terms in rrif/T and m&/T. Substituting <5i, 52 from 
Eqs. (2.21), we arrive at Eq. (1.6) . In the next Section we confirm the result (3.5) by 
evaluating the two-point correlation function of the stress-energy tensor in the sound 
mode channel and identifying the pole corresponding to the sound wave propagation. 
In addition to confirming Eq. (3.5), this will allow us to compute the sound wave 
attenuation constant and thus the bulk viscosity. 

4 Sound attenuation in J\f = 2* plasma 
4.1 Correlation functions from supergravity 

We calculate the poles of the two-point function of the stress-energy tensor of the 
Af = 2* theory from gravity following the general scheme outlined in [32]. Up to a cer- 
tain index structure, the generic thermal two-point function of stress-energy tensor is 
determined by five scalar functions. In the hydrodynamic approximation, one of these 
functions contains a pole at uj = u(q) given by the dispersion relation (1.1) and corre- 
sponding to the sound wave propagation in Af = 2* plasma. On the gravity side, the 
five functions characterizing the correlator correspond to five gauge-invariant combina- 
tions of the fluctuations of the gravitational background. The functions are determined 
by the ratios of the connection coefficients of ODEs satisfied by the gauge-invariant 
variables. Moreover, if one is interested in poles rather than the full correlators, it is 
sufficient to compute the quasinormal spectrum of the corresponding gauge-invariant 
fluctuation. This approach is illustrated in [32] by taking Af = 4 SYM as an example. 
For a conformal theory such as Af = 4 SYM, the number of independent functions 
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determining the correlator (and thus the number of independent gauge-invariant vari- 
ables on the gravity side of the duality) is three. In the M = 2* case, the situation is 
technically more complicated, since we need to take into account fluctuations of the 
two background scalars. These matter fluctuations do not affect the "scalar" and the 
"shear" channels 8 , entering only the sound channel. In the sound channel, this will lead 
to a system of three coupled ODEs for three gauge-invariant variables mixing gravita- 
tional and scalar fluctuations. The lowest or "hydro dynamic" quasinormal frequency 9 
in the spectrum of the mode corresponding to the sound wave gives the dispersion 
relation (1.1) from which the attenuation constant and thus the bulk viscosity can be 
read off. 

In this Section, we derive the equations for the gauge- invariant variables for a generic 
finite-temperature Pilch- Warner flow, i.e. without making any simplifying assumptions 
about the parameters of the flow. Then we solve those equations in the regime m^/T <C 
l,mj/T< 1, where the background is known explicitly. Solving the equations involves 
numerical integration. The final result is given by Eqs. (1.6), (1.7). 

4.2 Fluctuations of the non-extremal Pilch- Warner geometry 

Consider fluctuations of the background geometry 



where (more precisely, cf G , cf G ), cxbg, Xbg are the solutions of the equations of 
motion (2.8), (2.9). To simplify notations, in the following we use ci, c 2 to denote the 
background values of these fields, omitting the label " BG" . 
For convenience, we partially fix the gauge by requiring 



This gauge-fixing is not essential, since we switch to gauge-invariant variables shortly, 
but it makes the equations at the intermediate stage less cumbersome. We orient the 
coordinate system in such a way that the X3 axis is directed along the spatial momen- 
tum, and assume that all the fluctuations depend only on t,xs,r. The dependence of 
8 See [32] for classification of fluctuations. 

9 Gapless frequency with the property lim ui(q) =0 required by hydrodynamics. 



a 



®BG + , 

Xbg + $ , 



(4.1) 



X 




(4.2) 
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all variables on time and on the spatial coordinate is of the form oc e~ luJt+tqX3 , so the 
only non-trivial dependence is on the radial coordinate r. 

The fluctuations can be classified according to their transformation properties with 
respect to the 0(2) rotational symmetry in the x\ — X2 plane [8], [32]. The set of 
fluctuations corresponding to the sound wave mode consists of 



htti h aa — h xlXl -\- h X2X2 , hi X3) h X3X3 i 0, ip . 



(4.3) 



Due to the 0(2) symmetry all other components of h^ v can be consistently set to zero 
to linear order. It will be convenient to introduce new variables H a , H tz , H aa , H zz by 
rescaling 



hu — c\ H tt , ht z — c\ B. tz , h aa — c\ H aa , h zz — c\ H z 



(4.4) 



We also use Ha = H aa +H zz . Expanding Eqs. (2.8), (2.9) to linear order in fluctuations, 
we obtain the coupled system of second-order ODEs 



Hi + H' tt {Incjcl)' - H' u (lnd)' - \ U°\ H tt + u 2 H u + 2coq H tz 

C l V C 2 



8 (&P_ ^ + dV () 

3 \da dx 



H' t ' z + H' tz (\n^) +\cuqH aa = : 
\ Ci / c 2 

Ka + Ka (mc lC 5)' + (H' zz - H' tt ) (\nclY + - 2 L 2 - q 2 ^) H a 
3 \da dx J 



K z + H' zz (\n Cl 4y + (K a -H' tt ) (mc 2 )' + l 



+ q 2 -t(H tt -H a 



+ 8 (&P_ 

3 \da dx 



lu 2 H zz + 2uq H tz 
= 0, 



0" + 0' (In Cl cl)' + \ol BG (H u - Hu)' + \U 2 - 



i (d 2 v , d 2 v , 

oaox 



6 V da 2 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



(4.9) 
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r + i>' (In + \ X 'bg (Hu - H tt )' + 1 (c 2 - g 2 4 

_ 1 / <9 2 P 



1> 



+ 



<9 2 P 



(4.10) 



^ =0, 



2 <9x 2 

where the prime denotes the derivative with respect to r, and all the derivatives of the 
potential V are evaluated in the background geometry. In addition, there are three 
first-order equations 



ir Ull- 
Cl 



Hi, 



+ q 



C 2 



H' tz + 2[\n^\ H 



Cl 



Hz 



+ 8w (3a' BG + x'bg *l>) = , 



+ -f w #L -q H aa - 8q (3a' BG <p + x'bg = 



(lncic 2 )'^-(lnc3)'^ + 



w 2 + 2o;g fZ" ta + q 



2 u l 



2 v-^ft — H aa ) 



+ 4 



9a 



dV \ ,,, 

+ <9x V ~ 8 ^ 3a ^ G ^ + X ' BG ^ = ' 



(4.11) 
(4.12) 

(4.13) 



4.3 Gauge-invariant variables 

A convenient way to deal with the fluctuation equations is to introduce gauge-invariant 
variables [32]. (Such an approach has long been used in cosmology [41], [42] and in 
studying black hole fluctuations [43].) Under the infinitesimal diffeomorphisms 

x i*^ x n + £i* (4.14) 

the metric and the scalar field fluctuations transform as 

^^-VW(a, (4.15) 
i> -> i) - V x X bg f a , 

where the covariant derivatives are computed in the background metric. One finds the 
following linear combinations of fluctuations which are invariant under the diffeomor- 
phisms (4.15) 

Q 



Zh =4— H tz + 2 H zz - H a 



1 - 



q 2 c'jCi 
ui 2 c' 2 c 2 



L) 2 & 2 



a 



BG 



(lnc|) 

Xbg 



H n 



4%, ±J -aa i 



(4.16) 



(lnc|)' 
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Using Eqs. (4.5)-(4.13), one finds the new variables Z H , Z^, Z^ satisfy the following 
system of coupled equations 




(4.17a) 
(4.17b) 
(4.17c) 



where the coefficients depend on the background values C\, c 2 , olbg, Xbg- (The coeffi- 
cients are given explicitly in Appendixes B, C, D.) Eqs. (4.17) describe fluctuations of 
the background for arbitrary values of the deformation parameters rrib/T, m,f/T. 

The analysis of Eqs. (4.17) is simplified by switching to the new radial coordinate 
(2.10). Asymptotic behavior of the solutions to Eqs. (4.17) near the horizon, x — > + , 
corresponds to waves incoming to the horizon and outgoing from it, i.e. for each of 
the gauge-invariant variables we have Zh, Z^, Z^ oc x ±luJ . We are interested in the 
lowest quasinormal frequency of the "sound wave" variable Z H . To ensure that this 
frequency is indeed the hydrodynamic dispersion relation (1.1) appearing as the pole 
in the retarded two-point function of the stress-energy tensor of the Af = 2* SYM, one 
has to impose the following boundary conditions [32] 

■ the incoming wave boundary condition on all fields at the horizon: Zh, Z^, Z^ oc x~ lu) 
as x — ^ 0_)_ j 

■ Dirichlet condition on Z H at the boundary x = 0: Z H (0) = 0. 

The incoming boundary condition on physical modes implies that 



Z H (x)=x- iu Z H (x), Z+{x)=x- iu Z 4 ,{x), Z i ,(x)=x- i "Z f (x), (4.18) 



where Zh, Z^, Z^ are are regular functions at the horizon. Without the loss of gener- 
ality the integration constant can be fixed as 



Then the dispersion relation (1.1) is determined by the Dirichlet condition at the 
boundary [32] 



Following [44], [8], a solution to Eqs. (4.17) can in principle be found in the hydrody- 
namic approximation as a series in small u, q (more precisely, u/T <^ 1, q/T <^ 1), 



1 . 



(4.19) 



Z H (x) 



= 0. 



(4.20) 



X-t-l- 
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provided the background values ci, c 2 , olbg-, Xbg are known explicitly. Here we con- 
sider the high-temperature limit (2.14) discussed in Section 2.1, and expand all fields 
in series for Si <C 1, 5 2 <C 1. Accordingly, we introduce 



Z ° + 5? Z? + 5 2 2 Z 2 °) + *g + 51 Z{ + 5 2 2 , 



(4.21) 



where the upper index refers to either the leading, oc q°, or to the next-to-leading, 
oc q 1 order in the hydrodynamic approximation, and the lower index keeps track of the 
bosonic, Si, or fermionic, <5 2 , mass deformation parameter. Eqs. (4.18), (4.21) represent 
a perturbative solution of Eqs. (4.17) to first order in uj, q, and to leading nontrivial 
order in Si, <5 2 . From (4.19), the boundary conditions at the horizon are 



7° 



x->0 + 



1, z? 
0, z\ 



o, 



z 2 ° 



0, z\ 



x^0 + 



0, 

0. 



(4.22) 



The Dirichlet condition at the boundary (4.20) becomes 



7° 



X^l- 



X^l- 



= 0, z\ 
= 0, z\ 



X-t-1. 



X->1_ 



0, z 2 ° 
0, z, 1 



X-t-1- 



0, 

0. 



(4.23) 



X->1- 



We also find it convenient to parametrize the frequency as 



f2 



/3 2 <'2 



2g 



1 + (3{ Si + /3 2 r & 



(4.24) 



In the absence of mass deformation (Si — S 2 — 0), Eq. (4.24) reduces to the sound 
wave dispersion relation for the j\f = 4 SYM plasma [10]. The parametrization (4.24) 
reflects our expectations that the conformal j\f = 4 SYM dispersion relation will be 
modified by corrections proportional to the mass deformation parameters Si, S 2 - Our 
goal is to determine the coefficients (3%, f3[, /3 2 by requiring that the perturbative 
solution (4.21) should satisfy the boundary conditions (4.22), (4.23). 
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Using the high-temperature non-extremal Pilch-Warner background (2.16), param- 
eterizations (4.21) and (4.24), and rewriting Eqs. (4.17) in the radial coordinate (2.10), 
we obtain three sets of ODEs describing, correspondingly: 

■ the pure Af = 4 physical sound wave mode (Si = 0, S 2 = 0); 

■ corrections to pure Af = 4 physical sound wave mode due to the bosonic mass de- 
formation (Si 7^ 0, S 2 = 0) ; 

■ corrections to pure Af = 4 physical sound wave mode due to fermionic mass defor- 
mation (Si = 0, S 2 7^ 0). 

In the remaining part of this subsection we derive equations corresponding to each 
of these three sets. 

4-3.1 Sound wave quasinormal mode for Af = 4 SYM 

Setting Si = S 2 = in Eqs. (4.17), (4.21), (4.24) leads to the following equations 

x(x 2 + l)^r + (! - 3 * 2 ) 4^ + 4 ^o° = , (4.25) 
dx A dx 

rl 2 7 l J7 1 

x(x 2 + l) 2 -~ty — (x 2 + l)(3x 2 - 1) + 4a;(x 2 + 1) Z\ 

2 dZ» 4 ( 426 ) 

The general solution of Eq. (4.25) is 

Z ° = Ci(l - x 2 ) + C 2 [(x 2 - 1) lux - 2] , (4.27) 

where Ci, C 2 are integration constants. The condition of regularity at the horizon and 
the boundary condition (4.22) lead to 

Zl = 1 - x 2 . (4.28) 

Notice that the boundary condition (4.23) is automatically satisfied, as it should be, 
since (4.24) with Si = 0, S 2 = is the correct quasinormal frequency of the gravitational 
fluctuation in the sound wave channel for the background dual to pure Af = 4 SYM 
plasma. Given the solution (4.28), the general solution to Eq. (4.26) reads 

Z] = C 3 (1- x 2 ) + C 4 [(x 2 - 1) In x - 2] , (4.29) 

where C3, C4 are integration constants. Imposing regularity at the horizon and the 
boundary condition (4.22) gives 

Z] = . (4.30) 
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Again, the boundary condition (4.23) at x = is automatically satisfied as a result of 
the parametrization (4.24). 



4-3.2 Bosonic mass deformation of the Af = 4 sound wave mode 

Turning on the bosonic mass deformation parameter 5\ (while keeping 5 2 = 0) and 
using the zeroth-order solutions (4.28), (4.30), we find from Eqs. (4.17), (4.21), (4.24) 



3x 3 (x 2 - 1)' 



— x(x 2 — 1) a-y — , 



rl 2 7° d 7° 

x 2 (x 4 _ 1} _ x{x 2 _ 1)(3a , 2 _ i) 5£i + 4 x 2 (x 2 - 1) Z° 
dx 2 ax 



+ 192a; 2 



4x(l-x 2 ) " ! 



dx 



+ (1 +X 2 ) ft! 



+ 16x(x 2 - I) 3 (432) 



+ 32(x 4 -l)(x 2 -l)' 



+ 8x 2 (x 2 - 1) = 0, 



3x 3 (x 2 - if 



dx 2 



+ 3x 2 (x 2 -l) 2 
dZ« 



dZl 
dx 



- 2V3x 2 (x 2 - l) 2 — * - v^x 2 - 1) 
dx 



+ 3x 3 Z\ 



2(x 2 - 1) 



dai 
dx 



x ati 



0, 



(4.33) 



d 2 Z\ 



x 2 (x 4 - 1)(1 + x 2 
- 192x 2 (l + x 2 ) 



dZ\ 



- x(x 4 - l)(3x 2 -1)^L + 4x 2 (x 4 - 1) Z\ 
dx 1 dx 



Ax ( x 2 - 1) _ ( x 2 + 1) ai 



2 



dx 



dx 



x 2 (x 2 - l) 2 Z\ 



Z\ 



^v^x 2 



dcti 



2x(x 2 - l)(3x 2 + 1) — - - (1 + x 2 ) 2 ai 
dx 



32 

7! v 



dai 



-^=(l+x 2 ) 2 (x 2 -l) 3 (^) -i|x(x 2 + 3)(x 2 -l) 3CMl 



dx 



dx 



JL {x 2 + 3)(x 2 _ 1)x 2 ^ _ _^ x 2 (x 4 _ 1) ^ = , 



(4.34) 



where functions Ai(x), aii(x) are given by Eqs. (2.17), (2.19). 
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4-3.3 Fermionic mass deformation of the N = 4 sound wave mode 



Similarly, turning on the fermionic mass deformation parameter 5 2 and leaving 5i — 

we get 



12x 3 (x 2 



' dx 2 



+ 12x 2 (x 2 - 1) 



3x 2 (x 4 - 1) 



3x(x 2 -l) X 2 = 0, 
d 2 Z° 



dx 



dZ° 



+ 9x z Z° + 8(x 2 - l) 2 dX2 



dx 



(4.35) 



— f - 3x(x 2 - l)(3a; 2 - 1) ^1 + 12x 2 (x 2 - 1) Z\ 
dx 2 dx 2 



l 6x ( x 2_i) 3(1 + x 2 ) X 2 



-48a; 2 
+ 32(a; 4 -l)(a; 2 -l) 2 

^2 ■» d2Z ' 



dx2 



7° 



(4.36) 



dx J 



\2x\x 2 - iy 



+ A8x(x 2 - if 



dZ}, 



dA 2 
dx 



+ 24x 2 (x 2 - 1) (% = 0, 



Sv^a; 2 ^ 2 - if 



3a; 2 (x 4 - l)(l + x 2 ) 



dx 

d 2 Z\ 
dx 2 



V3(x 2 - 1) 



dZ 1 



(4.37) 



= 0, 



3a;(a; 4 - l)(3x 2 - 1) =?- + 12x 2 {x 4 - 1) Z\ 



-48x 2 (l +x 2 ) 



16x(a: 2 -l)^-3(l + a: 2 )x2 
dZ° 



- 2V3x(x 2 — l) 3 -j- 2 - + 4v / 3V(:r 2 - l) 2 Z 



- 32v / 3V 



8x(x 2 -l)(3a; 2 + l 



32 

7! 



dx2 

dx 
2 



4 



3X2(1 + x 2 ) 2 



(4.38) 



7° 



(i + I .)V-i)' (f) -Wa^W-D 

- Sv^a; 2 + 3) (a; 2 - l)a; 2 (% - 8v / 3a; 2 (a; 4 - 1) /3 2 r = , 
where functions A 2 (x), X2(^) are given by Eqs. (2.18), (2.19). 



dA 2 
dx 



5 Solving the fluctuation equations 

In this Section, we provide some details on solving the boundary value problems for 
the bosonic and fermionic mass deformations of the M = 4 SYM sound wave mode, 
discussing in particular the numerical techniques involved. We start with solving 
Eqs. (4.35)-(4.38) subject to the boundary conditions (4.22)-(4.23). 
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5.1 Speed of sound and attenuation constant to O in Af = 2* plasma 

Here we solve Eqs. (4.35)-(4.38) subject to the boundary conditions (4.22)-(4.23). 
Notice that the coefficient fi\ can be determined by imposing the boundary condition 
(4.23) on the perturbation mode Z%. The coefficient 0% can then De extracted by 
solving for the mode Z\ subject to the boundary condition (4.23). 

For the purposes of numerical analysis it will be convenient to redefine the radial 
coordinate by introducing 

Near the horizon (x = 0) we have y = x 2 + 0(x 4 ), while the boundary (x = 1) is 
pushed to y — > +oo. The computation proceeds in four steps: 

■ First, we solve Eq. (4.35). Applying the arguments of [32] to Eq. (4.35), we find that 
the appropriate boundary condition on Z° at y — > oo is 

Z\~y~^ as y^oo. (5.2) 

Eq. (5.2), along with the requirement of regularity at the horizon, uniquely determines 
the solution Z^(x). 

■ Second, we solve Eq. (4.36). The solution is an analytic expression involving integrals 
of the solution Z^(x) constructed in Step 1. Again, the regularity at the horizon plus 
the horizon boundary condition (4.22) uniquely determine Z%(x). The coefficient f3% is 
evaluated numerically after imposing the boundary condition (4.23). 

■ Third, we solve Eq. (4.37). The boundary condition 

4~y~ 3/4 as y^oo (5.3) 
uniquely determines the solution. 

■ Finally, we solve Eq. (4.38). The regularity at the horizon and the horizon boundary 
condition (4.22) uniquely determine Z\[x). Then the coefficient 0^ is determined 
numerically by imposing the boundary condition (4.23). 

Having outlined the four-step approach, we now provide more details on each of 
the steps involved. 

5.1.1 Step 1 

Asymptotic behavior near the boundary of the general solution to Eq. (4.35) regular 
at the horizon is given by 

Zl = A\y~^ + ■■■ + B\y~^ + • • • , (5.4) 
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where A^, £>° are the connection coefficients of the ODE. Rescaling the dependent 
variable as 

Z°^(l + y)-^ g^y), (5.5) 
we find that the new function g^{y) satisfies the following differential equation 

d 2 9j, y + 2 dg^ 9 3 /7 7 y \ _ 

dy 2 2y(l + y) dy 16y(l + y) 2 512y(l + y) 2 2 1 1,4' 4' ' 1 + y) 

(5.6) 

Imposing the regularity condition at the horizon, one constructs a power series solution 
near y = 

»* = <« + + 4) » + + sik) * 2 + 0(!/3) • (5 - 7) 

The integration constant g° is fixed by requiring that, as Eq. (5.2) suggests, 

g^~0(l) as y^oo. (5.8) 

For the numerical analysis, we had constructed the power series solution (5.7) to order 
0(y 19 ). Then, for a given we numerically integrated from some small initial value 
y = y in , using the constructed power series solution to set the initial values of g^(yi n ), 
gipiVin)- (This is necessary as the differential equation has a singularity at y — 0.) We 
verified that the final numerical result is insensitive to the choice of yi n as long as yi n 
is sufficiently small (we used y in = 10~ 15 ). Then we applied a "shooting" method to 
determine g®. The "shooting" method is convenient in view of the asymptotic behavior 
(5.4) near the boundary: unless g° is fine-tuned appropriately, for large values of y, 
9ip{y) would diverge, g^ oc y x l 2 as y — > oo. We thus find 

gl « -0.02083333(4) , (5.9) 

where the brackets indicate an error in the corresponding digit. We conjecture that 
the exact result is 

9i = ~^- (5-10) 
(Note that a formal analytical solution to Eq. (5.6) is available. It allows to express 
g° in terms of a certain definite integral. Numerical evaluation of the integral confirms 
the result (5.10).) 
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5.1.2 Step 2 



The solution to Eq. (4.36) regular at the horizon and satisfying the boundary condition 
(4.22) is 



Z°(x) = -(x 2 - 1) r z o(x) - - [(:r 2 - 1) Ins - 2] l^x) , 



where 



X z o (x) 



3x 2 (a; 2 + 3 + 21nx) 
4(1 + x 2 ) 2 



P V 2 + 



x , (z 2 -l)\nz-2 



6z 2 



-16z(z 2 -l)^ + 3(z 2 + l) X 2 



o z(z* - l)(z 2 + l) 2 



+ A(z 4 - l){z 2 - if 



dz ) 



+ 6z(z 2 - 1) ; 



dA 2 

dz 



l b z o(x)=- 



3x 2 



2(1 + x 2 ) 2 



P V 2 + 



f 

Jo 



dz 



z (z 2 + if 



X 



Qz< 



lQz(z 2 -l)^ + 3(z 2 + l) X2 



dxA' 



dz J 



+ 4(z 2 - l)(z 2 - if 
We explicitly verified that 

lim l z0 (x) = (9(1) 

a;— > 1 — 2 

Thus the boundary condition (4.23) becomes 



+ 6z(z 2 - if 



dA 2 

dz 



lim l 7 o(x) = 0. 



Numerically solving Eq. (5.15) for f3%, we find 

4 



&2 



3tt 



x 0.9999(5), 



(5.11) 



(5.12) 



(5.13) 



(5.14) 



(5.15) 



(5.16) 



where we factored out the value — 4/37T for the coefficient (3% obtained earlier from the 
equation of state (see Eq. (3.5)). 
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5.1.3 Step 3 

The general solution to Eq. (4.37) has the form 

(^/g \ ( I f x cos (-^arctanh z J 
— arctanh x ) x < C\ / dz — x 
2 J I QJo z 3 

Q 2/ 2 -,\ I o / 2 ,\ ^X2 „ 

8, (* - 1)^ + 8(3 -1)— -3, X2 



^/g \ |- 1 /"^ s i n ( -^arctanh z 



+ cos arctanh a; x < C 2 / ^ — x 

2 / | 6 ./ z d 



(5.17) 



8,V-D^ + 8(, 2 -l)%-3,X2 



where Ci, C 2 are integration constants. For generic values of these constants we have 
asymptotically 

Z\ ~ sin ( ^arctanh x\ {A%{\ + ■■■)+ B%{\ - xf' A + • • • ) 

/ /3 \ (5 - 18) 

+ cos ^-arctanh x \ {A% (1 + • • •) + - xf A + • • •) 

where B^, A c ^, B^ are the ODE connection coefficients. The integration constants 
Ci, C 2 should be chosen in such a way that the boundary conditions for the matter 
fields A^ = 0, A 1 ^ = are satisfied. 

To make numerical analysis more convenient, we introduce functions J- S {x), J- C {x) 

by 

Z^(x) = (1 — a; 2 ) 3 / 4 <jsin ^^^arctanh x^j ^(x) + cos ^^^arctanh x^j JF c (:r) 

(5.19) 

Redefining the radial coordinate as in Eq. (5.1) and using Eq. (4.37) we find that T s {x), 
J- C {x) satisfy the equations 



dJ- s 3 ^ 



s 



dy 4(1 + y) 

(5.20) 

3 (11 y \ 4, . dg^ 

128 2Fl U'4 ;3; TT^J + ^-3 (2/ + 1) ^ 

23 



cos ( ^arctanh y 



yV2(l + y)3/2 



dT c _ 3 
dy 4(1 + |/) c " 

(5.21) 

3 v ( 7 7 .o. y \ , A f , i\ ^ 



sin ( ^arctanh y 



2/1/2(1+^)3/2 



— 2 F 1 ( -, -; 3; + a. - -(y + 1) 

_128 V 4 4 ! + s vy ; 

Asymptotics of the solutions J- S {x), J- C {x) for y — > oo are 

* - 0(1) + ^ (l/ 3/4 ) , 
^ c - 0(1) + (y 3/4 ) . 

Accordingly, the initial conditions for the first-order ODEs (5.20) and (5.21) should be 
chosen in such a way that the coefficients of the leading asymptotics in (5.22) vanish. 
(This choice guarantees that the matter fluctuations of the mode do not change the 
fermionic mass parameter of the dual gauge theory.) Near the horizon, power series 
solutions to Eqs. (5.20), (5.21) are 

r, = f, + {-\g% - g) y 1 ' 2 + \f, y + Ofe 3 ' 2 ) , (5.23) 



128 4" c 8 " \ 32" L 8192 1536 



0(y 6 ) , 
(5.24) 

where /°, /° are integration constants, and g° is chosen as in Eq. (5.10). We use a 
"shooting" method to determine /°, /°: these initial values should be tuned to ensure 
that JF C , JF, remain finite in the limit y — > oo. We find 

/° ^0.01964015(5), 

( O..ZD J 

/ c ° w - 0.01743333(5) , 
where the brackets symbolize that there is an error in the corresponding digit. 

5.1.4 Step 4 

The solution to Eq. (4.38) regular at the horizon and obeying the boundary condition 
(4.22) reads 

Z\(x) = 2 ~^A 25x(x) + ^ 2jx(x) + ^ (a; 2 - 1) In* I°>(z) , (5.26) 

where the functions X^i(x), 2*i(x), 2T|i(^) are given explicitly in Appendix E. We 
verified that 

Km 2»i = (9(1) , Km 2|i = 0(1) , (5.27) 
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and thus the boundary condition (4.23) translates into the equation 



lim 2* i = . (5.28) 



To determine the coefficient 0%, we solve Eq. (5.28) numerically using the value of (31 
computed in (5.16). We find 

/3 2 r = (3j « 0.9672(1) . (5.29) 

5.2 Speed of sound and attenuation constant to O (5 2 ) in M = 2* plasma 

We now turn to solving Eqs. (4.31)-(4.34) subject to boundary conditions (4.22)-(4.23). 
The computation is essentially identical to the one in the previous Section, thus we 
only highlight the main steps. 

5.2.1 Step 1 

Using the radial coordinate defined by Eq. (5.1), we find that the asymptotic behavior 
near the boundary of the general solution to Eq. (4.31) is given by 

Z°, = Ay^ + --- + By^ log „ + ..., (5.30) 

where A^, £>° are the connection coefficients of the differential equation. Introducing 
a new function g^{y) by 

3js(l + v)- 1/S, *(v)> (5-31) 
we obtain an inhomogeneous differential equation for 

d 2 g lj> 1 dg<t, 1 



dy 2 + y dy 4y(l + y) 2 9<t> 96y( 



kw'Hv&TTi)- - (5 ' 32) 



Near the horizon one can construct a series solution parametrized by the integration 
constant g° of the solution to the homogeneous equation 

g * = 9 ° + (95 + i 9 *) 9 + (" m ~ s 9 *) 92 + ° (93) ■ (5 - 33) 

The integration constant g^ is fixed by requiring that 

g</,~0(l) as y^oo. (5.34) 
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Using the "shooting" method we obtain 



gl « -0.08333333(3) . 
We conjecture that the exact result is 



1 

12 



(5.35) 



(5.36) 



(Note that a formal analytical solution to Eq. (5.32) is available. It allows to express 
(?9 in terms of a certain definite integral. Numerical evaluation of the integral confirms 
the result (5.36).) 

5.2.2 Step 2 

The solution to Eq. (4.32) regular at the horizon and satisfying the boundary condition 
(4.22) has the form 

Zl{x) = 8(x 2 - 1) l^o(x) - 8 [(x 2 - 1) In x - 2] l b z o(x) , (5.37) 



where 



_ x 2 (3 + x 2 + 2lnx) f x (z 2 - l)lnz-2 

M w - - 4(1 + x 2 ) 2 " ^ + y dz ^(^_i)(i + Z 2)2 x 



2Az 2 



da i 



-Az(z 2 - 1) — - + (z 2 + 1) ai 



+ 4(z 4 -l)(^-l) 2 ('^ +6^ 2 -l) 3 ^ 2 



^?(^) = - n,i ~ 9V? Pi + I dz — - X 



dz f 



2(1 + a; 2 ) 



2Az 2 



Z (z 2 + 1)3 



-4z(z 2 - 1) — - + (z 2 + 1) ai 

dz 



dz 



+ 6z(z 2 - 1) J 



dAn 



+ 4(z* - l)(z 2 - V 2 < dX2 ^ ^ ' 



dz 



We have verified that 



tonics) = 0(1), 



and thus the boundary condition (4.23) translates into the equation 



(5.38) 



(5.39) 



(5.40) 



lim lL(x) = 0. 



(5.41) 
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Solving Eq. (5.15) numerically for fil we find 



32 

(3\ « x 1.00000(1), 

IT 2 



(5.42) 



where we factored out the value — 32/ir 2 obtained from thermodynamics (see Eq. (3.5)). 
5.2.3 Step 3 

Asymptotics of the general solution to Eq. (4.33) near the boundary y — > oo is given 
by Eq. (5.30) with different coefficients A\, B\. Rescaling the dependent variable 



Zl = {l + y)~ 1/2 G,(y), 
we find that the function G<p satisfies the following differential equation 

1 „ 



d 2 G ( f > 1 dG^ 
dy 2 y dy Ay(l + y) 2 
1 



G s + 



2^(1 + y) 2 



96y(l + y) 2 
0. 



2-^1 ( - ; t:! 3; 



2'2' + 



One can construct a series solution regular near y = 

^ o 1 o 



— G± + 



(5.43) 



(5.44) 



(5.45) 



The integration constant G° is fixed by requiring that 

G^ ~ 0(1) as y — > oo . 



(5.46) 



Using the "shooting" method and the value of the constant g° given by Eq. (5.36) we 
find 

Gj « 0.059(0) . (5.47) 

5.2.^ Step ^ 

The solution to Eq. (4.34) regular at the horizon and obeying the boundary condition 
(4.22) reads 



Z\{x) 



2(x 2 



3V2 



2(x 2 - l)lnrr 
3i72 



(5.48) 
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where the functions I z i(x), T zl (x), T c zl (x) are given explicitly in Appendix F. We 
checked that 

lim X zl ~ 0(1) , lim XU ~ 0(1) . (5.49) 

x^l_ 1 x— »1_ 1 

In view of Eq. (5.49), the condition (4.23) becomes 

lim 2* i = . (5.50) 



X-t-l- 



To obtain the coefficient /3f , we solve Eq. (5.50) numerically using the value of /3% given 
by Eq. (5.42). We find 

(3\ = Pi « 8.001(8) . (5.51) 

This completes out computation of the coefficients f3\, (3%, (3f, 0^: they are given, 
respectively, by Eqs. (5.42), (5.16), (5.51), (5.29). 



6 Conclusion 

In this paper, we considered the problem of computing the speed of sound and the 
bulk viscosity of Af = 2* super symmetric SU(N C ) gauge theory in the limit of large 
't Hooft coupling and large N c , using the approach of gauge theory/gravity duality. 
The computation can be done explicitly in the high temperature regime, i.e. at a 
temperature much larger than the mass scale m h and mj of the bosonic and fermionic 
components of the chiral multiplets, where the metric of the dual gravitational back- 
ground is known. Our results for the speed of sound and the bulk viscosity computed 
in that regime are summarized in Eqs. (1.6), (1.7), (1.8). It would be interesting to 
extend the computation to the full parameter space of the theory as well as to other 
theories with non-vanishing bulk viscosity. It would also be interesting to compare our 
results with a perturbative calculation of bulk viscosity in a finite-temperature gauge 
theory at weak coupling. 
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A Energy density and pressure in M = 2* gauge theory 

Energy density and pressure of Af = 2* SYM theory on the boundary dM$ of a 
manifold M$ with the metric (2.7) can be related to the renormalized stress-energy 
tensor one-point function as follows 

^ = M«VW, (A.l) 

P = v^iV s {T xixi )<f™, (A.2) 

where is the unit normal vector to a spacelike hypersurface E in CIM5, o is the 
determinant of the induced metric on E, and is the norm of the timelike Killing 
vector in the metric (2.7). The renormalized stress-energy tensor correlation functions 
are determined from the boundary gravitational action (with the appropriate countert- 
erms added) in the procedure known as the holographic renormalization. Holographic 
renormalization of M = 2* gauge theory on a constant curvature manifold was studied 
in [12]. Using the results for the renormalized stress-energy tensor one-point functions 
[12], one finds 

8 = 96^^ ( 18P ~ 9/3?1 " 12 *°* 10 + 36/3n/)l ° " + 36/5 ^ ' (A,3) 

' -e 4 « (6/3 + 9p 2 n + I2xlx w - 36p n Pio + 16*& " ^PuO , (A.4) 



967tG 5 

where the parameters /3, £, pio, pu, Xo, X10 are related to physical masses and the tem- 
perature for generic values of rrib/T and uif/T (see Section 6.4 of [31]). In the limit 
rrib/T <C 1 and m,f/T <C 1, matching the high temperature Pilch- Warner flow back- 
ground (2.16) with the general UV asymptotics of the supergravity fields we find [31] 



2, e« = 2 1 '\T 6*) 
\ n 2 7T J 



41n2 . A 8 . 

Pio = 01, pu = — di, (A.5) 

7T 7T 



\/2x 2[r(|)]- 

XO — — /QXl2 °2, XlO — 2 



[r(I)]- 



7T^ 



Using Eq. (A.5), from Eqs. (A. 3), (A.4) we obtain (to quadratic order in 5i, 5 2 ) the 
energy density and the pressure given in Eqs. (3.1). 
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B Coefficients of Eq. (4.17a) 

Ah(x) = 3uj 2 c 2 c\c\ ^—C2Ciq 2 c[ — 2c 2 c^q 2 + 3cj 2 c 2 c 2 ^ , 
B H (x) = u 2 c 2 c 1 [rjuj 2 c\c x cl - 42c?g 2 c' 2 3 + 9c\c x q 2 c 2 cl - %c\c\q 2 c 2 V 
+ %c\(? x q 2 d x V - 3c 2 c\q 2 d 2 2 d 1 + 9a; 2 c^c' 2 2 c'^j , 

C H {x) = w 2 ^-3w 2 4c^g 2 cici + 9u 4 c' 2 2 4 - 8q 2 c 4 d 2 2 c 2 2 V + 36q 2 c 3 1 d£d 1 c 2 

- 24q 2 44 + 16q 2 cld 2 d 1 (%P - %q 2 c\c? c A 2 V + 3c 2 c\q A d x c 2 + Qd 2 c\q A 

- \1q 2 c x clc\c 2 - lhu 2 c 2 2 clq 2 cl^ , 

D H (x) = lQq 2 c\^- Cl c 2 + d lC ^j ^2A Cl uj 2 (a BG )'Vcl + 36c 2 2 u 2 d 2 (a BG )' c 1 + 3c 2 2 c 1 u 2 ( 



C2q 2 c 2 1 ^c' 1 - 3Qc 2 c l d 2 u J 2 {a BG )' - 2±c 2 q 2 c\V{a BG )> - 2q 2 c\d 2 ^- 
oa ' oa 



Eh{x) = 16g 2 c 2 ^-cic' 2 + c[c2^ ^8cicj 2 x' bg 'Pc 2 + 12u 2 c 2 x'bg c 'i c I + 3c 2 ciu; 2 c 2 

- c 2 q 2 c 2 1 ^c' 1 - 8c 2 q 2 clV X ' BG ~ 12 Cl c 2 Vx' BG c 2 - 2q 2 c\d 2 ^j . 
The prime denotes the derivative with respect to r. 



,0P 

dx 
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C Coefficients of Eq. (4.17b) 

= 12c 2 c 2 c 2 ^— C2Ciq 2 c[ — 2c 2 c\q 2 + 3w 2 c 2 c 2 ^ , 

= 12cic' 2 c 2 ^-c 2 ciq 2 c 1 - 2c 2 c\q 2 + 3w 2 c 2 c 2 ^ ^3cic' 2 + c^c-^j , 

= -Qcluo 2 d 2 d iq 2 Cl - l§2cluo 2 c\{a' BG ) 2 V - 3c\u 2 d 2 2 c 2 ^ 

- 30c 2 2 LU 2 c' 2 q 2 c 2 + 6q A c 2 clc' lC2 + 2q 2 c\c\df^ - 48c^ 2 c 2 a^ G c^ 

+ q\\c\d 2 ^d x + 40q 2 4c' 2 a' BG 4^ + ^A^Bd^k + 12gV 2 2 c? 
+ 18c' 2 2 c 4 2 to 4 + l92q 2 c 4 c 4 2 (a' BG ) 2 V , 

- 16c 2 c 2 c 2 o;' BG — g 2 - 64c2C 2 a4 G g 2 XBG^ - 8c 2 C2C 2 gY BG ^ 

- %c 2 c x a BG —c x q 

( &P\ 

= -c\c\uj 2 \%^ bg Vc 2 + c>- , 

F<j, = cic 2 uj 2 ^8a' BG Vc 2 + ^- c ic' 2 + cic 2 ^ ■ 

The prime denotes the derivative with respect to r. 
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D Coefficients of Eq. (4.17c) 

„2 2 / / „ „ 2 / o '2 2 i o, ,2 / 2 



A,/, = 12c 2 c 1 c 2 ^— C2Cig c x — 2c 2 c 1 g + 3cj c 2 c 2 J , 
B^p = 12cic' 2 C2 ^— C2Ciq 2 c[ — 2c 2 c\q 2 + 3cj 2 c 2 c 2 ^ ^3ciC 2 + c'^^ , 

= 2^3q 2 cl4c' 2 ^c' 1 + 6q%clc' lC2 - 30c^ 2 c 2 2 g 2 c? - 64c^ 2 c 2 ( X ' BG ) 2 P 



- 48c^ 2 c 2 x^c 2 |- - 6 C ^ 2 c 2 cW 2 Cl + 6g 2 ctc 2 (c 2 ) 2 |^ + Uq 2 c\c\{x' BG ) 2 V 
+ 8? 2 c?x' BG c^c; + 40^4^^ - 9c^ 2 c 2 2 c 2 ^ + 18c 2 2 c^ 4 



+ \2q\d 2 fc\\ , 



= -2c 2 c 2 (-6c 2 2 c 2 J^g 2 + 9c' 2 2 c 2 2 -^u 2 - l92c 2 2 c\a' BG q 2 x' BG V 

dV dV 

+ \92c\ol BG u?x' BG V - 8c 2 CiXbgt^ c W 2 - 16c 2 c 2C 2 g 2 XB G ^ 

o / ^ / 2 , w ' 3 / 2 ^<-> / 2 / 2 

~ 2 ° 2Cl 9a9x Cig 72c 2 C 2 a BGT^ - 72c 2 c 2Ci«BGT^;<? 

o, / 3 2 / <9P\ 

+ 24d 2 4u 2 X ' BG —j , 

= -clcluo 2 h,—d 2 + 8c 2 x' BG V\ , 

= -c 2 cicj 2 ^3^c' 2 + 8c 2 Xs G P^ ^cic 2 - c^c^j . 
The prime denotes the derivative with respect to r. 
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E Structure functions of the solution (5.26) 



zLi(*) = 



3x 2 (3 + 21nx + x 2 



+ 



(l + x 2 ) 2 
4x 2 (7 + 7x 2 + 2x 4 



2x 2 (x 2 + 3f 



(l+x 2 ) 3 



lnx — ln(x 2 + 1) 



(l+x 2 ) 3 
2AV3z 2 (l + z 2 ) 



P v 2 + 



x , (z 2 -l)lnz-2 
dz^—, — — x 



o z(z±-l)(l + z 2 y 



-16z(z 2 -l)^ + 3(z 2 + l) X2 



J7° 

3z(z 2 -l) 3 -p + 6z 2 (z 2 -l) 2 Z\ 



+ 48z 2 [-8z(3z 2 + l)(z 2 - 1) + 3(1 + z 2 ) 2 X 2 

dz 

- 24z(z 2 + 3)(z 2 - l) 3 ^ - 16(1 + z 2 ) 2 (z 2 - l) 3 



Zl{z) 
dxi 



6x 2 



2x 2 (9 + 6x 2 + x 4 ) 



(l+x 2 ) 2 ^ 2 (l+x 2 ) 3 
24V3z 2 (l + z 2 ) 



+ / dz 



o ^(1 + ^ 2 )' 



dz 



-x 



16z(z 2 -l) ^ + 3(z 2 + l) X2 



dZ 9 ° 



-3z{z 2 -\f ^ + 6z 2 (z 2 -l) 2 ^2 



+ 48z 2 



-8z(3z 2 + l)(* 2 -l)^ + 3(l + z 2 ) 2 X2 

az 



4 



7° 



dAo 



-24z(z 2 + 3)(z 2 -l) 3 " 2 " 



- 16(1 + z 2 ) 2 (z 2 - 1) ; 



dz 

_ _ 6x 2 r _ 2x 2 (9 + 6x 2 + x 4 ) 

X ZlV X J — , 9^2 ^2 / 1 , „9\3 ^2 



d%5 



dz 



(l + x 2 ) 5 
24>/3^ 2 (l + z 2 



(l + x 2 ) J 



/ dz— 

Jo z(l + 



z 2 y 

4 



-16z(z 2 -l)^ + 3(z 2 + l) X 2 



dz 



_ 3 ^ 2 _1) 3 ^l 2 +6 ^ 2 (z 2 -l) 2 Z 2 ° 



+ 48z 2 



dX2 



-8z(3z 2 + 1)(^ 2 - 1) ^ + 3(1 + z 2 ) 2 X 2 

dz 

- 24z(z 2 + 3)(z 2 - l) 3 ^ - 16(1 + z 2 )\z 2 - 1) 



7° 

( dx2 



X 



dz 



dz 
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Structure functions of the solution (5.48) 

x 2 (3 + 21nx + x 2 ) r / 

^ZjW- (1 + X 2)2 ^ + [ 

4a; 2 (7 + 7a; 2 + 2x 4 ) 



2x 2 (x 2 + 3) 2 . l w 2 , ^ 
3(1+ a; 2 ) 3 ^-3^ + 1) 



3(1 + x 2 f 
96VSz 2 (l + z 2 ) 



PI 



x , (z 2 -l)lnz-2 
dz-—. — — x 



o z(z 4 -l)(l + z 2 ) 3 
Z\ 



-4z(z 2 -l)^l + (z 2 + 1) ai 

az 



^2 _ 1} 3 *£l + 2 ^2 _ 1)2 Z 



+ 192z' 



dai 



-2z(3z 2 + l)(z 2 - 1) — - + (1 + z 2 ) 2 ai 

az 



-8z(z 2 + 3)(z 2 -l) 3 ^ " 1 ■ - 2 - ? '- 2 ' 



dz 



-16(l + z 2 ) 2 (z 2 -l) J 



dai 

dz 



2*j(x) 



2a; 2 



r 2a; 2 (9 + 6a; 2 + a; 4 ) 



(1+a; 2 ) 2 ' 1 3(1 + a; 2 ) 3 

96v / 3z 2 (l + ^ 2 ) 



dz- 



4z(z 2 - 1) ^ + (z 2 + 1) «i 

az 



o z(l + z 2 Y 
Zl 



X 



- z(z 2 - 1) ; 



dZ\ 

dz 



2z 2 (z 2 - l) 2 Z\ 
da i 



+ 1922' 



-2z(3z 2 + l)(z 2 - 1) -± + (1 + z 2 ) 2 



-8z(z 2 + 3)(z 2 -l) 3 ^ ■ - 2 ^<~ 2 - ; 



dz 



-16(l + z 2 ) 2 (z 2 -l) J 



T c z} {x)=- 



2x 2 



(1 + a; 2 ) 2 



r _ 2a; 2 (9 + 6a; 2 + x 4 ) 



1- 

Jo 



dz 



7° 

dcti 
dz 

1 



96v / 3z 2 (l + 



3(1 + x 2 f 

, 4z ( z 2 _ i) ^1 + (^ + i) ai 

az 



z(l + z 2 ) 4 



_ ^2 _ 1)3 ^1 + 2z 2 ( ^2 _ 1)2 Z 



dz 



+ 192z z 



da i 



-2z(3z 2 + l)(z 2 - 1) ^ + (1 + z 2 ) 2 a, 



- 8z(z 2 + 3)(z 2 - l) 3 ^ - 16(1 + z 2 ) 2 (z 2 - l) 3 

dz 



d«i 

dz 
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